BOUNDARY LAYER FLOW PAST SYMMETRICAL 
AIRFOIL AT HIGH SUBSONIC SPEEDS 


By 

V. K. DOGRA 



A- 

m 


PF:i/977/ri 



DEPARTMENT OF AERONAUTICAL ENGINEERING 



INDIAN INSTITUTE OF TECHNOLOGY KANPUR 




JANUARY, 1977 




BOUNDARY LAYER FLOW PAST SYMMETRICAL 
AIRFOIL AT HIGH SUBSONIC SPEEDS 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

V. K. DOGRA 



r* 


to Che 


DEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JANUARY, 1977 











51TH. 




^ '' ^atP 1311 



CEETII’ICATE 





't'i 

ana!.'..®'/ 



0?his is to certify that the work "BOUNDARY 
MYEE EIOW PAST SYaHETEICAI AIEEOILS AT HIGH SUBSONIC 
SPEEDS" has been carried out under nQr supervision and 
has not been submitted elsewhere for a degree. 


Ail 

January 1977 (U .1,; Arora) 

Assistant Professor 
Department of Aeronautical Engineering 
Indian Institute of Technology 
KANPUE. 


i i" ■■ r ‘ !)U I B 

* i' - ■ ' ' d 

*'■■■ i • ■ ‘ ..f • 

‘ ^ c. -J.) , 

f iu ;t u> t 1 . . ■ i 

rcgulatio.K’ ■ , i lian 

institute oi i 

Dated. /, . 



ACKFOWIBDGEMEFTS 


Pirst and foremost, deep and sincere gratitude 
goes to Dr. If.L. Arora, for several reasons, among them 
being excellent and invaluable guidance, continued 
encouragement through the long tortpus work , and unfailing 
optimism. 


I thank Mr. R.N. Srivastava for his neat and 
accurate typing. 

finally my unbounded thanks to my friends who 
helped me a lot at different stages of the work and 
during my stay at I.I.T. Kanpur. 


January, 1977. 


V.K. DOGRA 



(i) 


TABLE OF GOI^TEMS 

SYNOPSIS (iii) 

NOIENCLATUEE (iv) 

LIST OP PIGUEES (vii) 

CHAPTER 1 INTRO LUC TION 1 

1 .1 General 1 

1 . 2 Literature Survey 2 

1 .3 Present Work 5 

CHAPTER 2 PROBLEM PORIiULATION 7 

2.1 Boundary Layer Characteristics 8 

2 .2 Thin Airfoils in Boundary layer Plow 1 1 

CHAPTER 3 INTEGRAL EQUATION APPROACH . 18 

3.1 Introduction 18 

3.2 Pormulation of Integral Equation 19 

3.3 Integral Equation for Perturbation 

Velocity Potential 25 

3.4 Integral Equation for Perturbation 

Velocity for Non-lifting Case 26 

3 .5 Simplification of Integral Equation 

for Perturbation Velocity 29 

3.6 Perturbation Velocity on the Surface 

of the Symmetrical Airfoil 31 



(ii) 


Page 


CHA.PTER 4 

IIETHOD OE SOLUTION 

34 

4.1 

.N-umerical Technique of Solution 

34 

4.2 

Method of Iteration 

36 

4.3 

Method of Parametric Differentiation 

38 

CHAPTER 5 

COMPUTATIONAL DETAILS, RESULTS AND 
DISCUSSION 

41 

5.1 

Introduction 

41 

5.2 

Method of Iteration 

41 

5.3 

Method of Parametric Differentiation 

42 

5.4 

Results and Discussion 

43 

5.5 

Conclusions 

44 

REPEREEOES 


46 

APPENDIX A 


48 


APPEEDIX B 


APPEEDIX C 



(iii) 


SYITOPSIS 


In tho present work the syimnetrical airfoil at 
zero incidonce in high subsonic bouncla,r7 layer flow has been 
considered. An airfoil in this kind of flow disturbs the 
boundary layer of flow field. The effect of perturbed 
boundary layer on the airfoil is incorporated by modifying 
its surface boundary conditions. 

lor high subsonic flows, the goTerning partial 
differential equation in terms of ' .. perturbation 

velocity potential, which is non-linear, is transformed into 
an integral equation with the help of G-reen's theorem. 

Using functional relationship for velocity distribution the 
double integral is reduced to an integral valid on the 
surface of the airfoil. The integral equation is replaced 
by a set of non-linear algebraic equations by quadrature. 

The set of non-linear algebraic equations is solved by 
simple iteration. In the second method solution to this 
set of algebraic equa.tions is obtained by applying tho method 
of parametric differentiation. This results in first order 
non-linear differential equations. The set of ordinary 
differential equations is solved numerically by .Hammings 
Predictor-Corrector Method. The pressure distribution is 
obtained for some typical symmetrical airfoils at zero 


incidence . 
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CHAPTER 1 


IHTRODUOTIOH 

1.1 Q-EHERAl 

The problem of calculating the pressure distribution 
on a wing of given shape is of central importance in aero- 
nautics. Generally the problem has been tackled considering 
the boundary layer on the airfoil or with the assumption of 
potential flow model. In the present work we calculate the 
pressure distribution on the airfoil tin steady compressible 
boundary layer flow, rather than a potential flow model. 

More specifically, instead of assuming that the initial flow 
velocity is constant, we specify that it varies with the 
co-ordinate normal to wing surface, and increasing to a 
constant (free stream) value within a finite distance <5, from 
the wing. This type of initial flow is called boundary layer 
flow. The velocity profile within this region of boundaiy 
layer flow is chosen to resemble that observed in viscous 
turbulent boundary layer. 

It is well-known that boundary layer flows are 
relatively more important for control surfaces than for the 
primary lifting surfaces. This type of flow is also created 
by some other agents, e.g., atmospheric boundary layer and 
propeller slip stream. 



The pressure distribution on the surface of an 
airfoil in inviscid flow is determined by the airfoil 
geometry, i.e., the thickness distribution, the shape of the 
camber line and the incidence. In the case of symmetrical 
airfoil at zero incidence, the pressure distribution is 
dependent only on the thickness distribution. But an airfoil 
in the boundary layer flow disturbs the stream-lines close 
to the surface due to the perturbation of initial flow 
boundary layer. This is equivalent to a change in the 
thickness distribution. Thus airfoil surface is modified 
by adding perturbed boundary layer thickness or some kind 
of average perturbed boundary layer thickness, e.g., 
perturbed displacement thickness. This analysis becomes 
more complicated for super-critical flows, because of shock- 
wave boundary layer interaction on the surface of the 
airfoil . 

1 .2 IiITEBATUBE SURVEY 

The shear flow past airfoil and wavy shaped wall 
has been considered by several authors. Tsion (17) consi- 
dered incompressible shear flow on symmetrical Joukowsky 
airfoils. There, velocity distribution in shear layer far 
ahead from airfoil is assumed some kind of linear function. 
The linear governing equation in terms of stre.am function 
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is generated and solved for perturbation velocities with the 
help of usual boundary conditions of the airfoil flow, 
lighthill (8) formulated steady compressible problem, deri- 
ving the inviscid perturbation equations and an asymptotic 
boundary condition. The appropriate boundary conditions are 
those of inviscid flow but with a mean boundary layer Mach 
number profile that goes to nonzero value at the wall. The 
incompressible problem for a stationary or travelling wavy 
wall has been treated by Benjamin (3) using full viscous 
perturbation equations. The effect of viscosity are included 
asymptotically near the wall. Inger and T/illiams (7) 
studied two dimensional steady compressible turbulent 
boundary layer flow past a slightljr wavy wall, in subsonic 
and supersonic cases , both analytically and experimentally 
in the Mach number range 0.8^ M 1 .8 at unit Reynolds 
niimbers of the order of 10° per inch. The analysis is based 
on a parallel shear flow model of undisturbed boundary layer 
and linearized treatment of the small perturbation field, 
where the solution is decomposed into slowly varying inviscid 
part, which determines the pressure and a rapidly varying 
viscous part that is important near the wall. The measured 
and predicted wall pressure distributions are in good 
agreement. Yates (19) developed linearized integral theory 
of the viscous compressible flow past a wavy wall, in that 
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two integr'als theories are derived for the calculation of 
pressure on a perturbed plane wall in a uniform main stream 
and in the presence of a viscous boundary layer. The 
integral theory may appropriately be called integral methods 
in that some kind of o-veraging over the thickness of boundary 
layer is employed j The first simple integral theory is 
derived from Von Karman momentum integral equation (19) while 
the second rational integral theory is derived by iteration 
from the asymptotics perturbation theory* Both theories yield 
the same functional form for the wp.ll pressure amplitude and 
phase angles in terms of a single similarity parameter, 
which is proportional to boundary layer thickness. The 
rational integral theory correlates better with experimental 
results than the simple theory. In recent years Ventres 
(18) presented a lifting surface theory based on a parallel 
shear flow model for steady, incompressible flows. The 
theory is intended to account approximately for the presence 
of a boundary layer. The method of Courier transfoms is 
used to calculate the pressure on a surface of infinite 
extent and arbitrary contour. Immediately above the surface 
is the region of sheared flow (the boundary layer), out- 
side of which the flow velocity is constant . The Fourier 
transform of the pressure on this surface is use^ to derive 
the shear flow equivalent to the kernel function of 
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classical potential flow lifting surface theory. The kernel 
function provides an integral relation between upv/ash at a 
given point on the surface and pressure everywhere on tho 
surface. It is found that shear lawyer decreases the lift 
curve slope; however, tho shear layer (whose thickness is 
constant along tho wing cord) has little effect on the centre 
of pressure. 

Another approach to the wavy wall problem is to 
divide the shear layer into discrete increments of uniform 
potential flow. The problem of each layer is solved exactly 
and solutions for a,djacont layers are matched, Tho limiting 
case of two layers was treated by Anderson and Fung (2). 
Zeyd-el (21 ) developed the theory for an arbitra.ry number of 
layers but did not present numerical results. 

1 .3 PEBSEFT WOBK 

The present work is an attempt to extend tho 
integral method of Yates (19) for boundary layer flow past 
wavy wall at supersonic speeds to boundary layer flow past 
two dimensional a,irfoils at high subsonic speeds. In 
Chapter 2 the perturbed displacement thickness is evaluated 
by perturbing momentum integral equation of compressible 
boundary layer about mean boundary layer flow i.d. parallel 
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shear flow of the undisturbed boundary layer flow in 
terms of perturbation velocity. Then the airfoil surface 
is modified by adding the displacement thickness to it. 

In Chapter 3? following Oswatitsch-Gullstrand-Sproiter-Alksone 
approach the differential equation is transformed into a.n 
integral equation by the application of Green's theorem, for 
high-subsonic flow. Using functional relationship for 
velocity distribution, the double integral is reduced to on 
integral equation valid on the surface of the airfoil. In 
Chapter 4, the integral is converted to a system of non- 
linear algebraic equations. The solution to the integral 
equarion is obtained by iteration scheme as well as by 
applying the method of parametric differentiation to the 
set of non-linear algebraic equations. Finally, results 
and discussion are given in Chapter 5. Some typical symme- 
trical profiles at zero incidence have been considered and 
pressure distribution is obtained a.t high subsonic speeds. 
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CHAPTER 2 


PROBIEM EOEIvIUIATIOH 


The method of treating airfoils in compressible 
bo-undary layer flows, is presented here. It may appropriately 
be called integral method in that some kind of averaging over 
thickness of boundary layer is employed. The analysis is 
based on a parallel shear flow model of undisturbed boundary 
layer, which is called mean boundary layer and perturbation 
boundary layer flow. 

The approach is directly based on Von Karman 
momentum integral equation and follows along the lines of 
Von Karman-Pohlhausen technique (19) for calculating average 
boundary layer properties on airfoils, wind tunnel walls, 
etc . This type of integral approach has been developed and 
extended successfully by lees (9) to more complicated flows 
that include, for example shock-wave boundary layer inter- 
action. Also an attempt has been made by Yates (19) to use 
the integral theory for calculation of pressure on a wavy 
wall, in uniform main stream and in the presence of shear 
layer at steady supersonic speeds. 
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2.1 BQ-PITOAEY lAYER CHABACTBRISTIGS 

We assimiG steady compressible boundary layer flow 
of a perfect gas with unit Prandtl number^ slibjectGd to 
steady perturbations sufficiently small. Considering X-ojcis 
along the chord of the airfoil and Y-axis is directed normal 
to the chord as in Figure 1 . Then momentum integral equation 
for compressible flows (13) can be \\Titten as; 


G 



(2 + H - M^) 


f 


T W 
e G 


( 2 . 1 ) 


where U is velocity of fluid in the direction of X-axis, 
is density and ji is coefficient of viscosity in the 
boundary layerj and M is Mach number at the edge of the 
boundary layer. The subscripts e, w denote the quantities 
at the edge of the boundary layer and on the surface of the 
airfoil, respectively. The other parameters are dofinod as; 




II 

(1 - ^)ax 

0 




e = [ 

0 

?U f. U 

'^0 8 C 

6* 



( 2 . 2 ) 




3u 


w 
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where 6 is the displacement thickness, 6 the momentum 
thickness, H shape factor and wall stress. Also, from 
the momentum equation at the edge of the boundary layer, 
we have 


dU, 


XT, 


e dX 


1 dp 
dX 


(2.3) 


The main idea is to perturb (2.1) about the basic 
mean boundary layer flow. The meaa boundary layer is plane 
parallel shear flow in the X-diroction with uniform pressure 
p , and arbitra,ry latereil variations of density f*^(y), 
velocity U^(Y) and Mach niimbor M^(Y). To do so, wo introduce 
two assumptions. First, we neglect the direct perturbation 
of the wall stress in (2.1). Second, we assume that the 
shape factor, H, is constant for the perturbed as well as 
the mean boundary layer flow. Thus 


and 




6 




o 


H(eQ + ©^) 





(2.4) 


where the subscripts o, 1 denote mean and perturbation 
quantities respectively. It follows from (2.4) that 

= I (2.5) 

In the absence of perturbations j (2.1) becomes 
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dO 

o 

dX 




e 


Subtra-cting (2.6) from (2.1), we get 


J, 

H 


d 6* 9 dU 
1 _o c 

dX ^ dX 


(H + 2 



0 


( 2 . 6 ) 


(2.7) 



This demonstrate simple proportionality between the pertur- 
bation displacement thickness and pressure coofficientj a 
key result of the integral theory. 

* ^ 

The moan displacement thickness 6^ and mean 
momentum thickness 0^ are assumed equivalent to flat plate 
values at zero incidence and at fixed Reynolds numbers 
based on the chord of the airfoil. These arc evaluated 
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accordin^^ to method given by Stratford and Beavers (14) on 

A 

the flat -plate. In this method the ;^th power law for 
velocit,/ profile in turbulent compressible boundary layer 
is used for computing the values of 6* and . The following 

set of results are obtained for free-stream Reynolds number 

f, 

R^, of the order of 10 

6* = 0.046(1 + 0.8M^ c 

m2 (2.10) 

= 0.036(1 + c r^^/5 

and K = '^q/©q 

where is the free stream Llach number and c is chord of 
the airfoil. Bor free-stream Reynolds number R_> of the 

7 

order of 10 the following are the results 

6* = 0.028(1 + 0.8M^)®*^^ c R^^^^ 

/ (2.11) 
= 0.022(1 + 7 ^ c R^^/^ 

and H = 

2.2 THIN AIRFOILS IR BOURBARY LATER FLOW 


A two dimensional cartesian co-ordinate system 
is chosen with the origin at the wing leading edge with 
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the X-axis in the free stream direction and Y-axis normal 
to free stream> as shown in Figure 1. The co-ordinates X 
and Y are scaled with respect to airfoil chord. A non- 
dimensional perturbation velocity potential,^, may be 
defined as 

= lx = ^ > ^ <2-12) 

where U and V are dimensionless perturbation velocities with 
respect to free stream velocity in X and Y dirootion 
respectively. 

The equation of motion governing inviscid, non- 
heat conducting potential flow of a perfect gas around thin 
airfoils at transonic speeds can be written as (1) 

[1 - (y + 1)^X"^XX ■^^YY ~ ° (2.13) 

The upper and lower surfaces of the airfoil may bo 
expressed as (Figure 1 ) 


= f^.(X) + f^(X) - ax 

(upper surface) 

= - f^(X) + f^(X) - OCX 

(lower siirface) 


(2.14) 

(2.14) 


where the body shape f is made dimensionless with respect 
to chord of the airfoil | f^(X) represents wing thickness 
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distribution, f (X) is equation for the camber surface andci. 

c 

is airfoil angle of incidence . 

The pressure coefficient C^CX, Y) can be evaluated 
by using Bernoulli's equation, thus 


Cp(X,Y) = 


Y 


r ^1 
^ L 


Y 


+ 


hi 

2 


M ^(1 - 

OO 


q^(X,Y))l ^ - 1] 


(2.15a) 


oo 

where y is the ratio of specific heats (y = 1.4) and, q(X,Y), 
is defined as 


q(X,Y) = 


(2.15b) 


A first order approximation to pressure coefficient consis- 
tent with (2,13) is given by 


C (X,Y) = - 2^^(X,Y) (2.16a) 

P 

which on the surface of the airfoil reduces to 

Cp(X, 0+) = - 2^j.(X, 0+) (2.16b) 

The complete problem for the perturbation field 
can now be expressed as follows. Bellowing second-order 
boundary layer theory (16), wo solve for the inviscid flow 
subject to the boundary condition, on the airfoil surface 
plus the displacement thickness, i.e, 

Y = f + ^*/c (2.17) 



14 


whore c is chord of the airfoil. Thus, we must solve (2.13) 
subject to the boundary condition in the plane of the air- 
foil (Z = 0+) which can bo expressed as 

4. df ^ 10^1 ro iRo^ 

0+) = (2.18a) 

using (2.9) together with (2.16b), we obtain from (2,18a) 

* 

|f(X, 0+) = + g ; + 2 - 1«2 (2.18b) 

In addition following conditions exe to be satisfied. 

All the disturbances produced by the airfoil must 
die out uniformly v/ith increasing distance away both from 
the airfoil and any downstream wake that may trail behind 
it i . 0 * , at inf inity 

<5 2; and<$y 0 (2.19) 

B’or lifting surfaces, it is necessary that^^ 
and % remain continuous and finite at the trailing edges 
i.o. Kutta condition is satisfied. 

It shovild -be noted that (2.13) describes correctly 
small perturbation in subsonic as well as supersonic flows, 
because then nonlinear term becomes negligible. The type 
of equation is dependent on the sign of the coefficient of 
as follows 
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[1 - ( Y + 

(50 OG 


’>•0 elliptic (subsonic) 

= 0 parabolic (sonic) 

CO hyperbolic (supersonic) 


Tho above coefficient is equivalent in the high subsonic 
flow -thoory to (1 - M ), whore M is the local Mach number. 
The critical value of U associated with M = 1 is thus 


1 


M" 


U 


ot. 


cr 


(Y + 1) M 


OC 


The lowest free stream Mach number corresponding to which 
the local Mach number becomes unity on the surface of the 
airfcil is called the critical Mach number. The flow with 
free stream Mach number less than critical value is called 
the subcritical flow, and the flow with free stream Mach 
number higher than critical value is called the super- 
critical flow. The supercritical flow in general is not 
shock free. Hence while treating supercritical flow, we 
have to take the effect of shock wave boundary layer inter- 
action, which makes the problem much more complicated than 
the subcritical. However, we have considered here only 
transonic viscous subcritical flows. 

The steady small disturbance supersonic potential 


equation of inviscid, non-conducting isontropic flow around 
two dimensional airfoil is obtained by neglecting the 
nonlinear toim in (2.15) > wo obtain 
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^ XX "^2^X1 ~ ^ (2.20) 

whoro A = - 1 

0 ® 

The equation (2.20) is hyperbolic and it can bo solved subject 
to the bounda.ry condition (2.18b) in the piano of the airfoil 
(Y = 0+) and in addition boundary condition at infinity (2.19) 
The general solution of hyperbolic equation (2.20) 
is given by 


5 = Pd -:Ay) (2.21) 

For the perturbation velocity components on the surf.ace of 
the airfoil, wo have 


^^(X, 0+) = 

^5y(X, 0+) = 


( 2 . 22 ) 


Using (2.22) and (2.16b) in the surface boundary condition 
of the airfoil (2.18b) 


^p(x, 0+) 

or 



0 


(H + 2 - M ) — ^ 


dC^(X,0) 


CO 


dX 


Cp(X,0) 


M. 

dX 



c 


(H + 2 



(f 

dC^(X,0-) 

dF 


or the upper surfac 
- Cp(X,0-) = - 


e) 

o M 
^ dX 


(2*23) 


(for the lower surface) 

Equation (2.23) is simple first-order differential equation 
for surface pressure coefficient. 



In the absence of any boundary layer, we 
(2.23) reduces to the usual inviscid result 

C = — ~ (inviscid) 

p ^dX 

The general solution for (2.23) can be v/ritten 




where 




M 

00 



= (H + 2? 



18 


CHAPTEH 3 

IITEGEAL EQUATION APPROACH 


3.1 INTEODUGPIOl? 

To find the pressure distribution around wing in 
high-subsonic boundary layer flowj we adopt the integral 
equation method « The integral equation approach for inviscid 
transonic flows has been developed by Ostwatitsch (12), 
Gullstrand (6), aiid Spreiter and Alksne (15)* Eecently this 
approach has been further worked out and refined by Norstrud 
(11) and Nixon and Hancock (10) . 

In this approach the application of Green's thoorem 
transforms the governing non-linear partial differential 
equation into an integral equation. That gives the pertur- 
bation velocity potential at a general point in the flow 
field,, in terms of a line integral of linear function of the 
velocity potential around the airfoil surface and a double 
(surface) integral of non-linear term. Por the perturbation 
velocity potential in the plane of the airfoil, this double 
integral is reduced to the simple line integral. 
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3.2 gOBMULATIOW OF IFTSGEAL EQUATIOI^ 


3.2.1 'Transformation 


It is convenient to introduce the following 
transformation 


M^(Y + 1) 
( 1 - ) 


(3.1) 


X = X 

7 = (1 - y 

CO 

Substituting (3.1) into (2.13)? (2.18b) and (2.19)? the 
problem of boimdary layer flow past airfoils at high sub- 
sonic speed is reduced to finding the solution of the 
equation 

^xx ^77 “ •^x^xx (3.2) 

together with the following boundary conditions. 

In the plane of the airfoil 

df 


C(x, 0±) = i s + 


where 


and 


Ai 


B. 


(Y + 1) 
(1 - 

CiG 


<5 

(H + 2 - M^)] 


(3.3) 




(1 - 

0<3 



20 


At infinity, ^ZL, 0 

y 


jZL, 0.. are finite and continuous at the trailing edge (Kutta 

y 

condition) i 

The pressure coefficient in the transformed plane 
is obtained from 


Cp(x, 0+) 


(T + 1) 


0 C^CjOi) 


3.2.2 Annlioation of G-reen's Theorem 


The equation (3.2) may be transformed into an 
integral equation by application of Green's theorem. 

If o(x,y) and^^(x,y) are two continuous functions 
with continuous first and second derivatives in somo domain 
S, bounded by curve G, then Green's theorem states 

'0'[^i’G(x,y) U'^cr(x,y) - ar(x,y) v jx(x,y)]dS 
S 




2? n 


d n 


dC 


(3.4) 


2 

where rr is Laplacian operator (~o +- — o ) » whore, n, 

is inward drawn normal to curve C around which the integration 
must be considered in an anticloclw/ise direction. 

In ( 3 . 4 ) let cj(x,y) be identified with pertur- 
bation velocity potential 0, so that (3.2) gives 
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= 

= ^4 (3-5) 

AgairitScXx ,y ) is chosen as the elomentary source solution of 
Laplace ' s oquation 

0 ( 3 . 6 a) 

i.e. SlMff f = ln(r) (3.6h) 


where r = [(x ~ ^)^ + (y - with ^ and ^ being the 

running co-ordinates in x and y directions respectively. 

ifoWvi^(x, yj is singular at the point P(x,y). 

Also to the present order of approximation, 0(^>^) and 
its derivatives can be discontinuous across the slit 
(0 ^ 1^1, r= Oi) on the ^-axis, where the boundary 
conditions are to be applied. Sinco^(x, 

0(x,y), and their derivatives must be continuous, throughout 
the domain S, the point P(x,y) and slit (0^^ 1, 

^ = 0+) must bo excluded from S, This domain is shown in 
Figure 2. 

On substituting (3.5) and (3.6) in (3j4), wo 


obtain 
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2 



J 


^ 9 









(3.7) 


whoro 0^ is the pa,rt of the boundary 0 surrounding the 
point P(x,y) and is talten to be small circle of radius r^. 

is that part of the boundary C surrounding the | -axis, 
and 0^ completes the boundary C, thus denoting the outer 
limiting boundary of S, taken to bo a large circle, centre 
at the origin, of radius E. Those boundaries C. , C 
together with the sense of integration, are shown in Figure 
2. The integral over surface S can be obtained directly 
without mentioning a.ny surface distribution, but for lino 
integral each boundary will bo considered separately. 


3.2.3 Integral Over Boundary 0 

(i) Integral o-ver boundary surrounding the singular 
point P ; 

The radius of circle surrounding the point P(x,y) 
is r^ and so 

60 = r 60 

P 

taking 0 to increase in clockwise direction. Since 
oTL (x, y| I ,^ ) = In(rp) 


? 
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3 _ b 

dn “ 3 Tp 


and 


(x, y; I ) 

dn 




= + 


Thus, tho integration over , in limiting operation a.a 
rp — ^ 0 is written as 

lim J [lll(r ) ac 

rp-^0 p CJ n an p 

2 It 

= - lim 0 r — . r dO 
r->0 o^ ^p P 

= - 2 TC0(x,y) (3.8) 


because all terms in line integral over vanish except 
those, involving derivatives of^Jb(x, )? Tp— ► 0. 

( ii ) Integral over boundary ; 

Por thin airfoil with thin boundary layer on its 
surface, we can write 


= + 


1 


d 

dn 


on tho plane of tho wing 
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Thus 




w 


= 1 Wx,! i y, 0) 


- { r, 0 ) ] a 


where 


^ 0(^,0) = ^(f,o+) - ^(1,0-) 


^ 0 (f,0) 


_ ^ 'MM _ Mho+l ^0(r,o- ) 


■3? 


(3.9) 


(iii) Integral over boundary Q 


lor integral over , the perturbation velocity 
potential 0C^)^ ) will be of the form 0(| ) = 0(9 + 

for a lifting airfoil (for non-lifting airfoil the 9 term 
should bo omitted), whore (R,9) or'e the polar co-ordinates 
on 0^. Thus the line integral over Gq^Is a constant for 
lifting airfoil. Then 




iMLiU 

d 11 




d n 



= Constant (for lifting case) 


= 0 (for non-lifting case) 


( 3 . 10 ) 


This constant is equal to the strength of circula.tion 
around the airfoil and exact value of this constant is 
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immaterial since it will disappear in the application of a 
differential operator. 


3.3 IFIECtHAL equation fOR PEETUBBAIIOIT VELOCITY POTEITTIAl 


Fow substituting ( 3 . 8 ), (3.9) and (3.10) in (3.7) 


wo get 


oo 


i2^(x,y) = [JL(x,| 5 7,0)21 j^^(r,0)- J^(x, I 5y,0)Ai2^(|'^0)] 


+ ^ Jufu(x,yi ^ ,^) ^(^(f>^))d.S + Constant (3.11) 

The surface integral in (3.11) can be defined in the same 
fashion as given by Spreiter and Alksno (15). Thus for 
field points in the upper half plane (y >• 0) 


// 


EdS 


^ x-e ^ ^ 

llm •/ f ( f pandr+ f ( J PaE)d5' 

e-^0 o -ea o x+e ^ 


o 

+ r ( r Ed| )d^ 

yj r\n J 


( 3 . 12 ) 


^00 -OO 


and for lower half plane (yc. O) 


(f EdS 
S 


.‘Soil 

} 


o CIO 


* J ( J PdE)a^ 

-CO X+e ^ 


The double integral in (3.11) can be integrated by parts 
with respect to|^ , so that 
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oo 


0(x,y) = ^ J^[v5L (x, I iy,0)Aj2L(|',0) — 

o ^ 6 " 


o 


^ -^J jJj ^) Constant 


2 n 


(3c13) 


wherGj^(x,yj | »|' ) 


in [(x -f)2 + (y 


3.4 IKTE&RAl EQUATIOE' gOE PEETUBBATION VELOCITY PQR KOIT- 
LIETIFG CASE 


In the symmetrical flow problem in the plane of the 
wing (y =■- 0±) 

0{x,O+) = 0(x,O-) 

so that /^0 = 0 (3.14) 

90(x.O+) _ 5 ^!(x.O-) 

d y ■ d y 


Thus we obtain from (3.13) 



(3.15) 

with u = 0 

Eow differentiating (3*15) with respect to x, we 
obtain for 0Ax,j) or u(x,y) 
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1 f 


= It J 


[(x -| ) + jI 


A0 (^jO)di 


+ 2 “ 2% ‘I '2 ^ 


IX - , ^ , i .. i am 

[(X - + (y -J^)2]2 ^ 


( 3 * 16 ) 


The term ^(x^y) arises from the differentiation of the 
limits eround the singular point P(x,y) in the surface 
integral* The boundary condition (3»3)) for symmetrical 
airfoil gives 


>0(x.O. 

dy 


n dx ax 


(3.17) 


Substituting ( 3 . 17 ) in (3.16), we get 


00 

A, r df. 


AW df^ (x 

w) = -i ^ .1 V/] 


^ ^ f 

n 4 


u C 


f tu-f 


%)2! y2] ^ 2 


2n 'if 2 


•(f,p 


r ^2^ 2 


[(x -f )‘ + (y 


d| 


(3.18) 


The equation (3*18) can also be written as 


« Uj^(x,y) - Ug(x,y) + ^(x,y) - I(x,y) 


u(x,y) 


( 3 . 19 ) 
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where 


U3^(x,y) 


cO 

h T "h (X -f ) 

” o ^ [(x -f )2 + y2] 




(3.20a) 


is essentially the linear perturbation velocity, 


%(x,y) 


B, 


Tl 


1 au(.k.,o.) 


■( x 


0 


ff v-s2 2-| ^ 

K [(x -^ ) + y ] 


(3.20b) 


is the boundary layer contribution, and 


I 
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(x - I )^ - (y - ^ 

[(x )^ + (y -^ )^] ^ 


dl d^ 


(3.20c) 


is the non-linear contribution. The integral u (x,y) can 

s 

be further simplified by integrating by parts, thus 




B. 


% 


u -r) 


00 


[(X -ir + y^] 


u 


(f,o) ! 

o 




[ (x - f )^ - y^] 

[(x + y^]' 




It is assumed that the perturbation flow ahead of and 

behind the airfoil has very sma,ll effect on the values of 

u (x,y), hence the first term in above integral can bo 
s 

neglected . Thus 


Ug(x,y) 


B. 


71 


oo 

J 


u{F,0) 


[(x - 1)^ - y^] 
[(x - ^) + y ] 


(3.20d) 
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How from (3.18), (3.19) and ( 3 . 20d ) , we get 


OCJ 


u(x, 7 ) = 


It 


a? 


(X -I) 


[{X -f )2 + y2] 


d^ + 


ci6 

B 


% 


j ^(8 0) -.f-(x,7) 


o 


[(x - I )2 + y2] 2 



Rx -I )2 - (y -.r )2 1 

[(X -5 + (y 2 



( 3 . 21 ) 


3.5 SBlPLIBICAflOH OB IFTBGBAL EQUAflOH BOR PERTUBBAJIOH 
VELOCITY 


In solving (3.21) the difficulty arises due to the 
kernel of double integral. Therefore before proceeding 
further the kernel of double integral and hence integral 
equation is simplified. Bellowing Spreiter and Alksne (15) 
a functional relationship 


u(x,y) = 


u(x,0i^) 




( 3 . 22 ) 


is assumed between the velocity on the airfoil surfa.ce and 
the velocity directly above the surface of the airfoil. The 
expression (3.22) satisfies the boundary condition of zero 
disturbance at +C30. The parameter b(x) is found by using 
the condition kf irrotationality and the boundary condition 
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(3«3)i Diff erontiating (3*22) with respect to y and setting 
y = 0+, we have 

2u(x,0+) 


b(x) = - 


ou 

!v=:l 


y=0+ 


From the irrotationality condition 


3 u ® 

5?!, 


y=0+ 


Uoi 


The boundary condition (3.3) gives 

+ A. ^ + B- ‘|-^(x,0+) 
' dx ' Cx*^ 


ar 

(5^ 


7^2 

li, 


y=o+ 


(3.23) 


(3.24) 


(3.25) 


Substituting (3*24) and (3.25) into (3.23) > we get 

-2u(x,0+) 

b(x) 


(Y+1)M^ 6* ^ (l-M^) ^u(x,0+) 

2 -^ + i7r(s+2-ii^i -- ^ -p -- ] 

o J (Y+DmIb 


The typical calculation shows that (H + 2 - M ) f*\^0{10”^). 




OC 


u 


This when multiplied with (which is of the some order of 
^2 Ox^ 

magnitude as becomes one order higher as compared with 

dx 

the first term in the denominator, and therefore ce.n bo 


neglected . Thus 


b(x):i=i - 


2u(x,0+)(1 - ]£)^/^ 


(Y + 1) 




a 

^ae> , 2 
dx 


(3.26) 
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3.6 PERTURB AT lOU VELOCITY ON THE SURE ACE OP THE SYMI3TRICAL 
AIRPOIL 


Substituting relation (3.26) into (3.21), we obtain 


, at 7 = 0 


u(x,o) = Tr/ ^ (x -I) irf 


al 


(x )' 


CO 


u 


:(x,o) - n ff C(^ -g . ^l 

^ -V ^ [(=^ -?) ] 


dS a^ 

^ At • ^ ^ 


+ ^(x,0) 


(x ) 


2 2 


oo 


-kid i^(f,o-) 

-00 ^ ^ rfx (1 -'I- 


[(x (1 


+ f 0) -lCx....,.-. l .}f. _d| ] ^ 


(3.27) 


Since the airfoil is symmetrical, we have 
u(|,0+) = -u(|',0-) 

Using (3.28) in (3.27)? we get 


(3.28) 
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OD oo 


(x -tr 
> 

2 r2 


2 ^ CO JO r, 

+ ^(x,0) - lit / [ / u^(|,0) — C ? -. - - ? .,. )— - iX ^ 

2 ^ -ij 0 ^ [U -r)2 




(1 +^)2 ^ 


The (3.29) is further simplified as 


A. p df 

■u.(z,0) = ~ J (3^ t 


dl 


B. 


+ 


u 


d| (x ) 
o6 

( 5 , 0 ) - r 


4 * 


71 


i" i^(|j0) 


df 


(x -f )' 


--OC ^ 


E(Z) d 


where E(Z) = E( i--^-^ -) 


OO 


i f (x -§ 

J r / \2 . 2 


[(x -I ^ [1 + 




Bow writing (3*30) as below 


u(x,0) = Uj^(x,0) + Ug(x,0) + ^(x,0) - iCxjO) 


where 


Ut.(x, 0 ) = 


I 


df. 




% J df (x -?e) 


(3.29) 


(3.30) 


(3.31) 


(3.32a) 


The linear perturbation velocity on the surface of the 
airfoil 
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and 


00 

1 


u(P,0) 


AL 


(x -C ) 


(viscous contribution) 


(3.32b) 


06 


l(x,0) = 1 ^2b ' ( ' gl ' ^ ' (non-linear contribution) 

~ OO s . . 

(3.32c) 

The linear perturbation velocity Uj^(x,0) in (3.31) 
is singular at the leading edge for round nosed airfoils. 
Therefore, for more accurate solution this is corrected to 
new modified value u^(x,0) (see Appendix A) as 




Uj^(x,0) = 


P. 


+ Ut(x,0) 


OO 


-if t 2, i 

r-i + (SQ — 1) ] ^ 


where 


k 


(y + 1 ) and = V 1 - M 


00 


c£) 


Thus, substituting this modified linear value in (3.31), we 
get 

2 

= Uj^iXfO) + Ug(x,0) + ^(x,0) - I(x,Q) 


u(x,0) 


(3.33) 
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CHAPTER 4 


l/ETHOL OE SOIUTIOH 


4.1 HUaERIGAL TEOHRIQUB QE SOLUTION 

The integral equation (3.33) of Chapter 3 is solved 
by quadrature , by approximating the range of integration 
with a finite umber IT of discrete intervals and letting the 
unknown function u = u(x,0) be represented by a mean value 
within each element, Eor points ahead of the leading edge 
(x <<C0) or points aft of trailing edge (x >1), the expression 
for u as given in (3-33) is going to be identically zero. 

Thus it is assumed that the influence of the points upstream 
and downstream of the profile on the value of u(x,0) can be 
neglected. Thus (3.33) can be represented by the following 
system of non-linear algebraic equations . 


Ui = 


-pH P H 


where c . . and 6 
11 11 


i=1 i=1 

■i - 1, 2 Nf 

are defined as (5) ' 





(4.1) 
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^ Vt TaT + 8(1+A^) - 8 ] 

12(1+A )^ 

+ 12A(A^-1) In |Aj - A(1+A^) [(1+A^) + 12]2 


12 ( 1 -. b 2)^12 ?Bj 


^o^lir [d+B^)^ - (1+B^)^+ 8(1+B^) - 8] 


+ 12B(B^-1) In iBi - B(1+B^) [(1+B^)^ + 12]| (4.2a) 


A = 




(4.2b) 


J 


h-W^ 


a, -^)- 


X. - f . + TT 

1 2 


(4.2c) 


where 


A 

(step size) 


Vi^riting (4.1) in following form 


- 2Ui + 25i,, - 2 L 


ij 2 


u. IT 

77^ + 2 5~ 6. .u . 

2 IJ 0 


i = 1 ,2, . . . ,IT 


(4.3) 
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How to execute (4*3) numerically, the value of u 

versus x is to be considered as divided up into a number of 

equally or unequally spaced strips j over which u is assumed 

to be constant in steps. If 1 . is taken to represent the 

«] 

width of the jth strip of the divisions into which x-axis is 
partitioned j represents the mid-point of this jth strip 
and u. stands for the ordinate corresponding to the height 

vl 

of this stripj then at this particular jth step c. . and 6 . . 

1 J 1 J 

are defined as in (4.2a) and (4.2b). 

Thus a set of non-linear algebraic equations (4*3) 
can be solved by different iteration schemes. The solutions 
to these equations can also be obtained by the use of arti- 
ficial viscosity in various numerical calculation methods 
and one of these methods which^uso here is the method of 
parametric differentiation. This approach is due to Horstrud 
(11) applied to nonlifting problem. 

First we shall discuss the method of iteration 
scheme for present nonlifting problem and then the 
parametric differentiaticn method for. the same. 


4,2 METHOD OF ITERATION 

The integral equation (3.33) for perturbation 
velocity ■ can bo solved by the method of iteration. To 
do so we arrange (3*33) as follows. 
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u(x,0) = u^(x,0) + ^(x,0) - It;(x, 0) + Ug(x 50 ) (4.4') 


where = (I - 

oo 

I _ ( .a! (l,0) 

J _ 2b E 


and 


ii 


(Z)d^ 


oo 

This can he further expressed in the form of a qua,dratic 
equation 

u^(x,0) - 4u(x,0) + 4^Uj^(x,0) - l^(x,0) + Ug(x,0)‘|= 0 

(4.5) 

Equation (4.5) has the following roots 


u 


(x,0) = 2[1 ± 1 1 - (u^(x,0) - It(x, 0) + u^(x,0))| ] (4.6) 


Since the flow is high subsonic for which u(x,0) <C2, wo 
retain only the negative sign before the radical in (4.6). 
Thus we have 


u 


(x,0) = 2 [1 - - (u^(x,0) - I^(x 50 ) + u^(x,0))^ ] (4.7) 


First of all, Uj^ is computed numerically. 

Then, equation (4.7) is operated for u(x,0) by neglecting 
integrals l^(x,0) and Ug(x,0). Now, with this computed 
value of u(x,0) the integral l^(x,0) is computed and 
substituted in (4*7), keeping integral u^(x,0) equal to 
zero to get new value of u(x,0). With this new value of 
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u(x,0) tho integral l^(x,0) is computed again and this 

computed value is substituted in (4^7)) keeping u^(x,0) 

equal to zero. The above procedure is repeo,tod till the 

value of u(x#0) converges. Then with this converged value 

of u(x,0) the integral u (x,0) is computed and put in ( 4 . 7 ) j 

s 

that gives the final value of u(x,0); 


4.3 METHOD Of PARAIVIETRIC DIT'T'EBSNTIATIOU 

The method of parametric differentiation is applied 
for solving the non-linear algebraic equation (4.3) for the 
unknovm u^ (20). To do so a parameter jCt is introduced as 
follows 


F^(Uj_(u) ) ^i) 



- r 1 2 

I O 




H e • -u . 
fr^ ID D 


+ ST 6- . 1 = 0 

3-* 


(4.8) 


i = 1 ,2 , . . . 

By differentiating (4.8) with respect to tho parameter il, 
we got a system of first order non-linear ordinary differ- 
ential equations. Considering u^ independent of |^, because 
Uy is solution of (4.8) at p. = 0 and unknown u. = u. (p), 
we obtain tho set of ordinary differential equations as 
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1 

(1 - [I Uj^) 


rl 2 




0=1 


2 


du . 



H du. 

+ T. + -=^)] 

j=1 ^3. 0 


(4.9) 


1- 

Tho unknown vector can bo determined from differential 
equa,tions as the solution to the defined Cauc]:: 57 ''s problem 
at jl =-1, with the initio-l vector solution 


u.=Ut. a.t 5=0 (4.10) 

1 Jj . ^ 

1 

The condition for a unique solution of system (4.9) in the 
range 0 ^ {I 1 is the nonvanishing of the functional 
determinant or Jacobian of F^[u^(il), p ] with respect to u^^ 
for fioced p. Jacobian matrix of the system is given by 


J (^-j > ^2 > ' ' • 


d(]F g^) 
^ (u-j j . • . »''^J^) 


(4.11) 


and (4.11) is positive definite under the condition that the 
elements on the principal diagonal are positive, i.o. , 


(1 - p u^) > 0 

Uj; < , i = 1 ,2, ... ,17 

For a certain value of p- ^ max [u^^ ] however Jracobian m.atrix 
becomes singular and the solution curve of Cauchy problem 
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can split into two or more solutions. Such a point p = fl* 

in the interval 0 ^ jl 1 for v/hich det J(u-|, ~ ^ 

be treated as bifurcation point and the associated solution 

Uj|^ = u^(p ) is designated as the bifurcation solution. 

After obtaining numerically a tmique and real 

solution Uj^ = u^(p) for (4.9) at value p ^ 1 y which 

represents within some given error bound, the best possible 

of 

approximation to the condition minimum bifurcation, the next 

A 

step in the integration of (4.9) will involve crossing of 
singular Jacobian matrix. Therefore a continuous dependence 
of the solution on the initial data is no longer possible . 
Then^ method of steepest descent is applied for the system 
of algebraic eqijation (4.8). Por continuous flows the 
initial values for the steepest descent method are those at 
which the Jacobian matrix becomes singular. Prom the 
bifurcation point the iteration by steepest descent method 
is carried upto fi = 1 to get final velocity distribution. 

In the present case, Hamming’s predictor-corrector 
method in conjunction with fourth order Eunge-Kutta starter 
has been employed to integrate the set of differential 
equations . 


can 



CHAPTER 5 


GOMPHTATIOHAL DETAILS, RESULTS Affi) PISCUSSIOH 

5.1 IHTRODUOTIOR 

The computational details for the methods used in 
solving the system of equations (4-3) are given in Chapter 4. 
Results are computed and compared for HACA 0012 airfoil at 
different free stream Mach numbers by the method of iteration 
and parametric differentiation. The loading edge correction 
slightly lowores the value of linear coefficient of pressure 
at the loading edge of the airfoil Eigure 5 ”. 


5.2 METHOD OE ITERATION 

fy 

The flow chart is shown in Eigure O . listing of 
the computer program is given in Appendix B. Definition of 
the variables used in the program is given in the listing. 
The slope of the airfoil is denoted by E. The subroutine 
KARMAH ca,lculates the linear velocity distribution DBAR. 
This lines,r velocity is modified by applying leading edge 
correction and used to find the first approximation of U. 
Then, the subroutine IIPMAN is called, which uses the above 
U as the initial value and calculates the inviscid non- 
linear effects by iterating each now value of U until 

, ■ >' 

l.i. ■ 
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convergonce is achiovod. The subroutine G-IIJERT uses the 
converged inviscid value of U and calculates the fincil values 
of U in boundary layer flow* 

5*3 METHOD OR PARiaiETRIG RIgFEEBMPIATIOR 

The listing of the developed program for the 
problem under donsideration is given in Appendix 0. The flow 
chart describing the main stops of the computer progrrn is 
shown in Figure 4 ^. The different variables used for compu- 
tation are defined in the program listing. The equation 
for lACA 0012 o.irfoil is denoted by Y(A) . The linear velocity 
distribution UBR, is cnlcula.ted separately and fed to the 
main program. 

For applying the Hamming's predictor-corrector 
method the initial values o,ro oalculcated using the Runge- 
Eutta method. A step size H is assumed and the initial 
conditions for the Runge-Kutta method are taken at IJ- =0, 

U = UBR. The subroutine FOTO' calculates the velocity 
deriva,tivG F. In the subroutine - FUNG, the subroutine 
MITAl is called which calculates the influence coefficient 
B. In the Eunge-Rutta method, wherever the velocity 
deriva,tivo is required, the subroutine FUNG is called. 

In the Hamming's predictor-corrector method, the solution 
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vector USP is first predicted. The solution vector USP is 
then modified using the corrector formula. A convergence 
test is applied after each step to ensure that the values 
are within the proscribed error bounds . After calculating 
the final velocity distribution USH", at different AtCEW, a 
check is applied for the Ja,cobian matrix. 

5.4 RESULTS mi) DISGUSSIOK 

Results are obtained for the UACA 0012 airfoil at 
Mach numbers of 0.4, 0.5, 0.6, 0.7 and 0.72. It is observed 
that the results obtained in boundary layer flow are higher 
than that for a non-linear inviscid flow. Eigures 8,9,10,11 
and Eigure 12, show the comparison of the results obtained 
for coefficient of pressure. 

The presence of a boundary layer on the surface of 
the airfoil produces a distortion of the streamlines close 
to the surface equivalent to a, change in the thickness 
distribution. Consequently, the perturbation velocity of 
the outer inviscid flow on the modified surface of the air- 
foil, increases. This expla.ins as to v/hy the coefficient of 
pressure is more negative for the bound a.ry layer flow as 
compared to that of the non-linear inviscid flow. It is 
also observed that the difference between the two results 
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ai*e more pronotincGcl at higher Mach numbers because the 
displacement thickness^ increo.ses with an increase in the 
perturbation velocity < 

The results obtained by tho method of pa.rametric 
differentiation end the method of itorn.tion Figure T 
agree except near the leading edge of the airfoil at low 
Mach numbers j This is because tho number of stops taken in 
the method of iteration is twice that taken in tho method of 
paramotric differentiation. The slope of tho airfoil, at 
tho loading edge, changes rapidly and hence tho difference 
in the number of steps taken produces an appreciable diver- 
gence in tho results. 

Tho calculations were performed on the computer 
system IBM 7044. The order of magnitude of tho computer time 
used for getting the final results for the methods of iter- 
ation and paramotric differentiation are given below J 

Method of iteration 3~5 minutes 

Method of paramotric differentiation 5-8 minutes 

5.5 COMCLUSIOUS 

To conclude, v\ro observe that tho presence of a 
boundary layer on the surface of tho airfoil increases 
the negative coefficient of pressure as compared to that 
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for tho non-linear inviscid flow. The method of iteration, 
for calculating the problem of boundary layer flow past on 
a.irfoil at a high subsonic speed, is better suited than 
the method of parametric differentiation as the computer 
program is simple to construct and the computational times 
involved are less. 


V 
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APPEiTOIX A 


SIECtULAEITY corrections at IBADIHG EDGE OE THE AIREOIL 


OiiG main disadvantage of the linear theory is that 
the solution ofUj^x,0) as given by (3i32a)| is singular at 
the loading edge; It implies that a more accurate solution 
to that given by linear theory should be used in the region 
of the loading odgoj preferably without any singular 
behaviour. 

Now consider the governing differential equation 

(2.13) . . 
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X “XX 


with boundary condition on the surface of the airfoil 
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where f represents the surface of the airfoil. Applying 
the transformation (3.1) to (A. 2), we get 
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Wow from (2.15b) together with the transformation (3»l)} 
we get at the surface of tho airfoil 

q(xsO+) = t^+UL(x,0+)] [1 + (p^) f (A. 4) 

The exact flow remains finite everywhere so that at the 
nose of the airfoil 


U-r (X 5C + ) 


XI 


k , df 

p2 - dx 


<50 


Then (A. 4) can bo written in the form 


q(x,0+) = 


ih * [1 - w s ] 


(A. 5 ) 


where both the numerator and the denominator tend to 
infinity at nose. 

If the standard linearised solution Uj^(x,0) given 
by ( 3 . 32 a) is to be used as an approximation to modified 
linear velocity Uj^(x,0+) and since this linearised solution 
already incorporates an infinity at the leading edge then 
an approximation to q(x,0+) is taken to be 


k 


2 + Uj^(x,0+) 


6^ df , 2 ^ 
^ ^k dx ^ 


q(x,0+) 


(A. 6) 
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It is noted that only the thickness terin is retaiiiod in 
(A<6), This is done for two reasons ; first it is only the 
thickness term which gives necessary infinity at loading 
edgej second it is necessary to preserve the same denominator 
for both the upper and lov/er surfaces to ensure that Kutta 
Condition a.t the troAling edge is satisfied, since '5il(x,0+) 
already satisfies Kutta condition. 

An approximation to the 5(x»0+) may be v/ritten as 


q(x,0+) 


B2 


+ Ut(x,0+) 


(A.T) 


where Uj^(k,0+) can be regarded as the linear solution 
modified for leading edge corrections. Combining (A, 6) and 
(A.T), we obtain 
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(A. 8) 
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